Model of two-fluid reconnection 
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INTRODUCTION 

Magnetic reconnection is a fundamental physical pro- 
cess of topological rearrangement of magnetic field lines 
in magnetized plasmas during which magnetic energy 
is converted into kinetic and thermal energy. Recon- 
nection is of particular importance in the solar atmo- 
sphere, the Earth's magnetosphere, and in laboratory 
plasmas [3, 0, 0, 0). In hot plasmas, because of the low 
electrical resistivity, magnetic reconnection due to resis- 
tive dissipation alone is very slow. As a result, a simple 
single-fluid MHD description of the plasma is generally 
believed to be insufficient for the theoretical explanation 
of fast reconnection events. Instead, a two-fluid MHD ap- 
proach has been frequently used in recent studies of fast 
reconnection [H, 0, 0,13] . Most these studies have been nu- 
merical and experimental, while an ultimate goal of con- 
struction a comprehensive theoretical model of two-fluid 
reconnection has not yet been achieved. In this Letter 
we present a model of two-fluid magnetic reconnection, 
which serves this goal. 



TWO-FLUID MHD EQUATIONS 

We use physical units in which the speed of light c and 
four times tt are replaced by unity, c — 1 and An = 1. To 
rewrite our equations in the Gaussian centimeter-gram- 
second (CGS) units, one needs to make the following 
substitutions: magnetic field B — > B/-v/4vr, electric field 
E — * cE/ v47T; electric current j — ► V^Lrj/c, electrical 
resistivity i] — > ijc 2 /Att, the proton charge e — > y/Awe/c. 

We consider an incompressible, non-relativistic and 
quasi-neutral plasma, composed of electrons and protons. 
Using standard notation, the equations of motion for the 
electrons and protons are [H, 0] 

nm e [<9 t u e + (u e V)u e ] = -VP e - ne(E + u e xB) - K,(l) 
nm p [d t u p + (u p V)u p ] = - VP p + ne(E + u p x B) + K, (2) 

where n is the (constant) number density, and the sub- 
scripts and superscripts e and p refer to electrons and 
protons respectively. Here, K is the resistive frictional 
force due to electron-proton collisions that can be ap- 
proximated as K = n 2 e 2 r](u e — u p ) — —nerjj, where rj 



is the electrical resistivity [2, [B| . For simplicity, we ne- 
glect proton-proton and electron-electron collisions and 
the corresponding viscous forces. We also introduce the 
electric current j = ne(u p — u e ) and the plasma velocity 
V = n(m p u p + m e u e )/ p, where p = n(m p + m e ) = const 
is the plasma density. Taking into account m e <C m p , we 
find u p = V + m e j/nem p and u e = V — j/ne. We sub- 
stitute these expressions and K = — nerjj into Eqs. ([1]) 
and |(3J). We also substitute the electric field E from 
Eq. ((I]) into Eq. ©. We obtain 



E 



tq - V x B + j x B/ne - [VP e - {d 2 Jd 2 p )V P p ]/ne 
+ d\ [da + (W)j + (jV)V - (l/ne)(jV)j] , (3) 
nm p [dtV + (W)V] = —VP + j x B ^(jV)j, (4) 

where P = P e +P p is the total pressure, d e — (irie/ne 2 ) 1 ' 2 
and d p — (rrip/ne 2 ) 1 ^ 2 are the electron and proton iner- 
tial lengths. Eq. (J3j> is the generalized Ohm's law describ- 
ing the motion of the electrons. Eq. (j4|) is the plasma mo- 
mentum equation describing the motion of the protons. 
We note that the electron inertia terms, proportional to 
d 2 , enter both Ohm's law and the momentum equation. 
Although these terms are essential for fast two-fluid re- 
connection (as we shall see presently), they have been 
frequently neglected in the momentum equation before. 
We also note that V • V = and V • j = for incompress- 
ible and non-relativistic plasmas. 

We consider two-fluid magnetic reconnection in the 
classical Sweet-Parker-Petschek geometry, shown in 
Fig. [TJ The reconnection layer is in the x-y plane with 
the x- and y-axes perpendicular to and along the layer 
respectively and all z derivatives are zero. The thickness 
of the reconnection current layer is 2(5, which is defined 
in terms of the out-of-plane current (j z ) profile across 
the layer It can be shown that 28 turns out to be 
also the thickness of the electron layer, where the elec- 
trons are decoupled from the field lines. The length of 
the electron (current) layer is defined as 2L. The pro- 
ton layer, where the protons are decoupled from the field 
lines, has thickness 2 A and length 2L ex t, which can be 
much larger than 26 and 2L respectively. The values 
of the reconnecting field outside the electron layer (at 
x w 5) and outside the proton layer (at x « A) are 
about the same, B y « B ext up to a factor of order 
unity. The out-of-plane field B z is assumed to have a 
quadrupole structure (see Fig. [T]) [2j, yj, |4|. Also, the re- 
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FIG. 1: Geometrical configuration of the reconnection layer. 

connection layer is assumed to have a point symmetry 
with respect to its geometric center O (see Fig. Q} and 
reflection symmetries with respect to the x- and y-axes. 
Thus, V x (±x,Ty) = ±V x (x,y), V y (±x,Ty) = TV y (x,y), 
V Z (±X,TV) = V z (x,y), B x (±x,Ty) = TB x (x,y), 
B y (±x,Ty) = ±B v {x,y), B z (±x,^y) = -B z {x,y), 
jx(±x,Ty) = ±3x(x,y), j y (±x,Ty) = Tj y (x,y) and 
j z (±x, +y) = j z {x, y). The derivations below extensively 
exploit these symmetries and are similar to @, 0| • 

SOLUTION FOR TWO-FLUID RECONNECTION 

We make the following assumptions for the reconnec- 
tion process. First, rj is constant and small. Second, the 
reconnection process is quasi-stationary, so that we can 
neglect time derivatives. This assumption is satisfied if 
there are no plasma instabilities in the reconnection layer, 
and the reconnection rate is sub-Alfvenic, E z <C VAB ext . 
Here Va = B ext / \jnm v is the Alfven velocity. Third, the 
pressure tensors are isotropic, so that the pressure terms 
in Eqs. |3|) and (j4|) are scalars. 

Using Ampere's law and neglecting the displacement 
current, we find the current components to be j x = d y B z , 
j y = —d x B z and j z = d x B y — d y B x . The z-component 
of the current at the central point O is 

jo = (jz)o = (d x B y - d y B x ) o ps B ext /S, (5) 

where we use the estimates (d y B x ) <§; (d x B y ) and 
(d x B y ) ps B ext /S at the point O. The last estimate 
follows directly from the definition of S. 

Equation (0| for the plasma (proton) acceleration 
along the reconnection layer gives 

nm p (VV)V y + d 2 e (jV) Jv = -d y P + Jz B x - j x B z . (6) 

The y derivative of this equation at the point O gives 

nm p {d y Vy) 2 (l + d 2 e f/d 2 p ) ps 2B 2 ext /L 2 + 3o (d y B x ) . (7) 

Here we introduce a useful dimensional parameter 

7 = (d xy B z ) /ne(d y V y ) , (8) 

which measures the relative strength of the Hall term 
(jxB) z /ne and the ideal MHD term (VxB) z inside the 
electron layer. In the derivation of Eq. we use the 
estimate (d yy P) w {d yy B 2 /2) ext » ~2B 2 xt /L 2 , which 



follows from the force balance condition for the slowly 
inflowing plasma across the layer Q. 

Faraday's law V x E = — <9tB for the x- and y- 
components of a quasi-stationary magnetic field in two 
dimensions gives d y E z = —d t B x as and d x E z — d t B y ss 
0. Thus, E z is approximately constant in space, and from 
the generalized Ohm's law ^ we obtain 

E z = r)3z - V x B y + V y B x + (j x B y - j y B x )/ne 

+ d 2 e [j X 3 X V Z + jydyVz + Vxd X j Z + Vydyj Z 

- (jxdxjz + 3ydyjz)/ne] w constant. (9) 

The reconnection rate is determined by the value of E z 
at the point O, namely E z = rjj . We estimate j a below. 

Taking the second derivatives of the z-component of 
Eq. ([I]) with respect to x and y at the point O, we 
find (d xx V z ) = -{e x f/m p )[{d x B y ) - d 2 e (d xx j z ) ] and 
(d vy V z ) = (ej/m p )[(d y B x ) + d 2 e (d yy j z ) }. Using these 
expressions, we calculate the second derivatives of Eq. © 
with respect to x and y at the point O and obtain 

i 1Jo /6 2 a 2(d y V y ) oJo [l + d 2 e /S 2 } 

X [1 + 7(1 -di^/dl)], (10) 

rjjjL 2 » 2(d y V y ) [(d y B x ) -d 2 eJo /L 2 } 

x[l + 7(l-^/ d 2 } ] 5 (n) 

where we use (d x V x ) a = —(d y V y ) , (d x B y ) w j a , 
(d xx j z ) « -jo/5 2 and (d yy j z ) w -jo/L 2 . 

In Eq. j9]), the electric field E z is balanced by the MHD 
and Hall terms outside the electron layer, where the elec- 
tron inertia terms are unimportant. Therefore, 

E z » -V x B v (l-jx/neV x ) w (d y V y ) S B ext (l+% (12) 
E z « V y B x (l- ]y /neV y ) pb {d y V y ) (d y B x ) L 2 {\+^) (13) 

at the points (x « S, y = 0) and (.t = 0, y ps L) re- 
spectively. Here we use the estimates j x w {d xy B z ) 5, 
jy ps -(d xy B z ) L, V x ps -(9yUy) (S, U y ps (9 y F y ) L, 
-Ba; ps (d y B x ) a L and i? y ps i3 ea; t. 

The ratio of Eqs. {TTJ and HJ) gives (d y B x ) ps 
(<5 2 j /L 2 )(l + 2df/5 2 ), while the ratio of Eqs. (0 
and US]) gives (d y B x ) a ps B ext S/L 2 ps 5 2 j a /L 2 , where 
we use Eq. ([5|). Comparing these two estimates, we find 
5 > d e ; therefore, j < B ext /d e and £^ z < rjB ext /d e 

Next, we use the z-component of Faraday's law: ps 
—dtB z — d x E y — d y E x . Taking the c^y derivative of this 
equation at the point O, and using Eq. ([3|) for £^ and E y , 
after tedious but straightforward derivations, we obtain 

ps -r) [(d xyxx B z ) + (d xyyy B z ) ] + (1 - d 2 e ^f/d p ) 
x[(d y Bx)o{dxxjz)o + (d x B y ) (d yy j z ) ]/ne 
« r,ne(d y Vy) j/6 2 - (1 - d 2 7 /4) 

x[(9 y B x ) j /5 2 +i2/L 2 ]/ne. (14) 

To derive the final expression, we use Eq. © and the 
estimates {d xyxx B z ) a ps -{d xy B z ) /S 2 > (d xyyy B z ) 0l 
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(d XX j Z ) ~ -jo/5 2 , {dyyj Z ) C 



Note that Eq. flU} results in 7 < d 2 p /d 2 e . 

We estimate the proton layer half-thickness A as fol- 
lows. Outside the electron layer the electron inertia and 
magnetic tension terms can be neglected in Eq. ([6]), and 
we have nm p (VV)F 9 ~ —d y P. Taking the y derivative 
of this equation at y = 0, we obtain nm p [V x (d xy V y ) + 
(d y V y ) 2 ] » -{d yy P) « 25^/L 2 . Here the term 
V x {d xy Vy) is about of the same size as the term (d y V y ) 2 . 
Therefore, we find that {d y V y ) ext ~ Va/L outside the 
electron layer (but inside the proton layer). Next, in the 
upstream region outside the proton layer ideal single- fluid 
MHD applies. As a result, at x « A and y — Eq. © 



-jo/L 2 , (d x B y ) j a . occurs when d p « L ex t/V~S [1, H, 0, H] ■ Therefore, our 



reduces to i? z = 77^0 
AB ext V A /L. Thus, 
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(0„V„) ea! iAB e 



Va/L, A » r)j L/V A B e 



(15) 



Now we solve equations ©, (H)-©, (UOj) ([T5j) for un- 
known quantities j , S, A, L, 7, {d y V y ) , (d y B x ) and 
(d xy B z ) . We neglect factors of order unity, and we treat 
the external field B ext and scale L ext as model parame- 
ters. Recall that parameter 7, given by Eq. ([5]), measures 
the relative strength of the Hall term and the ideal MHD 
term in the z-component of the Ohm's law. Depending 
on the value of 7, we have the following solutions and the 
corresponding reconnection regimes. 

Sweet-Parker reconnection. When 7 < 1, both 
the Hall current and electron inertia are negligible, and 
the electrons and protons flow together. In this case, we 
obtain the Sweet-Parker solution: j w \fS B ext / L ext , 
Ez = V3o ~ V A B ext /VS, S « A w L ext /VS, L rj 
L ext , 7 ~ Sd 2 /L\ xt , (dyVy)o w (d y Vy) ex t « V A /L exU 
(d y B x ) Ri B ext I L ext \fS , (d xy B z ) « SB ext d v /L\ xt , 
where 5 = VAL ext /ri is the Lundquist number. Con- 
dition 7 < 1 gives 5 < L<1 xt /d 2 . Therefore, Sweet-Parker 
reconnection occurs when d p is less than the Sweet- 
Parker layer thickness, d p < L ext /\/~S @, H, 0|- 

Hall reconnection. When 1 < 7 < d p /d e , the Hall 
current is important but the electron inertia is negligi- 
ble. In this case, the solution is j « B 2 xt /r\neL = 



JneL = {d p /L)V A B e 
A ! 



Sd p B ext / LL ext , 

S w ner]L/B ext = LL ext /Sd p , A » a p , 7 
Sd 2 /LL ext , (d y Vy) » (d y V y ) ext ~ Va/L, {d y B x ) « 
ne^/I, = B ext L ext /Sd p L and (p xy B z ) w B 2 xt /nei]L 2 = 
Sd p B ext / L 2 L ext . These results are in agreement with 



previous theoretical findings 0,0, EEi. Condition 1 
dp/d e translates into Sd e d p / L ex t < L < Sd 2 /L e 



< 



7 < 

for the electron layer length L. Unfortunately, the ex- 
act value of £ cannot be estimated from Ens. (I51)-(IT51) 
in the Hall reconnection regime. In studies [7|, [9j, [ll| L 
was essentially treated as a fixed parameter. Here, we 
take a different approach and make a conjecture that the 
Hall reconnection regime describes a transition from the 
slow Sweet-Parker reconnection to the fast collisionless 
reconnection (presented below). Numerical simulations 
and experiment have demonstrated that this transition 



conjecture implies that the Hall reconnection solution is 
S ~ L 2 xt /d 2 , L ex t > L > d e L ex t/ dp, j rj B ex tL ex t/d p L, 
E z Rj {dp/L)VAB BXt , 8 « d p L/L ext , A w d p , 7 « L ext /L, 
(d y V y ) w {d y V y ) ext ~ Va/L, (d y B x ) rj B ext d p /LL ext 
and (d xy B z ) a w B ext L ext /d p L 2 . As the electron layer 
length L decreases from its maximal value L rj £ext to 
its minimal value L sa d e L ext /d p , this solution changes 
from the slow Sweet-Parker solution to the fast collision- 
less reconnection solution that is presented next. 

Collisionless reconnection. When d p /d e < 7 < 
d 2 /d 2 , the electron inertia and the Hall current are im- 
portant inside the electron layer and the proton layer 
respectively. In this case, the solution is 

jo ~ B ext /d e , S w d e , (16) 

E z =rij rj riB ext /d e = (L ext / Sd e )V A B ext 

w (A/Z^Besi w (d p /L)V A B ext , (17) 

L « V A d e d p /r) = Sd e dp/ L ext , (18) 

dp/de < 7 < d 2 p /d 2 e , A^d p , (19) 

(9 y n)o « »?A#y = V A L ext /Sd% (20) 
{d y V y ) ext ~ v/dedp = V A L ext / Sdedp ps Va/L, (21) 

(d y B x ) « B^/Vf^ - B ext L 2 e Xt /S 2 ded 2 p, (22) 
(d xy B z ) a fts Bextv/VAdedp = B ex tL ex t/ Sd e d p . (23) 

Apart from the definition of the reconnecting field -B e:c t, 
Eqs. (fT6)) -([T8 |) essentially coincide with the results ob- 
tained in [8] for an electron MHD (EMHD) reconnec- 
tion model. Note that the value of 7 or, alternatively, 
the value of the proton acceleration rate (9yV^) at the 
point O cannot be determined exactly. This is because 
in the plasma motion equation ([6]), the magnetic tension 
and pressure forces are balanced by the electron iner- 
tia term d 2 (y\7)j y inside the electron layer. The proton 
inertia term nm p (V\7)V y can be of the same order or 
smaller, resulting in the upper limit (d y V y ) < Va/L. 
Thus, inside the electron layer the magnetic energy is 
converted into the kinetic energy of the electrons (and 
into Ohmic heat), while the proton kinetic energy can be 
much smaller. However, in the downstream region y > L, 
as the electrons gradually decelerate, their kinetic energy 
is converted into the proton kinetic energy. As a result, 
the eventual proton outflow velocity becomes w Va [23| • 
These results emphasize the critical role that electron in- 
ertia plays in the plasma momentum equation (|4|). 

The collisionless reconnection rate given by Eq. (|17|). 
although being proportional to resistivity [24J, is much 
faster than the Sweet-Parker rate E z Rj V A B ex t/ \fS as 
long as S < L\ xt /d\. The solution lfIB j) -p5 ]) is valid pro- 
vided L ex t/d e <C S < L 2 xt /d e d p , which is obtained from 
conditions E z <C VAB ex t and L < L ex t- Thus, both fast 
collisionless and slow Sweet-Parker reconnection regimes 
can exist simultaneously, as found in simulations [l3j |. 
These simulations also found a hysteresis for transition 
between slow and fast reconnection regimes. This implies 
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that the transition Hall regime, during which the elec- 
tron layer length L decreases, may occur at Lundquist 
numbers other than S ~ L? ext jd 2 p , depending on the past 
history. Unfortunately, our stationary model cannot de- 
scribe time-dependent transition processes. 

It is known that the single-fluid MHD reconnection be- 
comes much faster when resistivity r\ is anomalously en- 
hanced by current-driven plasma instabilities [3, 0,0 Hi- 
Eq. (fT7|) shows that resistivity enhancement can consid- 
erably increase the collisionless reconnection rate as well. 
This enhancement can occur after the electric current 
value (j ) jumps up during the transition from the Sweet- 
Parker to the collisionless regime at d p f=a L ext / y/S, and 
could be a physical mechanism of very fast reconnection. 



electric field E z is supported by the Hall term (j x B) z /ne 
in the downstream region y > L. As a result, there 
are Hall-MHD Petschek shocks attached to the ends of 
the electron layer (25[, as observed in numerical simu- 
lations 21|. However, in these simulations a spatially 



localized anomalous resistivity was prescribed, resulting 
in a short layer length, while here resistivity r\ is assumed 
to be constant. 
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was supported by the NSF Center for Magnetic Self- 
Organization (CMSO), NSF award #PHY-0821899. 



DISCUSSION 



Let us compare theoretical results (fT6|) -(|23 |) for col- 
lisionless reconnection with numerical simulations and 
experiment. The estimates A sa d p for the proton 
layer thickness, S w d e for the electron layer thick- 
(d xy B z ) 5L w B ext for the quadrupole 



ness, B z w 

field, and u y w -j y /ne w (d xy B z ) L/ne pb v eA 
Bext I s/nme for the electron outflow velocity agree with 
simulations 0, i, S, [H, [H, E3, El • The estimates A « d p 
and B z rs B ext also agree with experiment Q. However, 
the experimentally measured electron layer thickness is 
about eight times larger than the model and simulations 
predict [3, [2(| ■ Three-dimensional geometry effects and 
plasma instabilities make direct comparison and exact in- 
terpretation of the experimental results difficult 0, [2(| ■ 
Our theoretical results are qualitatively consistent with 
recent numerical findings of an inner electron dissipation 
layer and of electron outflow jets that extend into the 
proton layer 15 , 3, 17 , HH ■ The electron layer length 

1/2 

L cc d e ex m e decreases with the electron mass, as 
in simulations [18], but the scaling law observed in these 
simulations was slightly different, L oc m^ 8 . Length L w 
Sd e d p / L ex t is generally much larger than both 6 ss d e and 
A sa dp, consistent with simulations HI, 3, 17 1. How- 
ever, if resistivity 77 becomes anomalous and enhanced 
over the Spitzer value so much that S ~ L ext /d e , then 
L can theoretically become of order of d p , as in sim- 
ulations [H, [ijj]- Our theoretical results for the pro- 
ton velocity V agree with simulations [l7j], which found 
the proton outflow velocity to be significantly less than 
Va and also found acceleration of protons in the decel- 
erating electron jets. Unfortunately, detailed quantita- 
tive comparison of our results to the results of kinetic 
numerical simulations is hindered because the simula- 
tions do not explicitly specify resistivity 77. Also, in the 
simulations the electron pressure tensor anisotropy was 
found to play a critical role inside the electron layer and 
, while in this study an isotropic pressure is 
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jets 

assumed and the electrons are coupled to the field lines 
inside the electron outflow jets. Thus, in our model the 
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